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QJ . Abstract 

Q 

Tlj" I We explore the prospects of low-scale leptogenesis in a class of supersymmetric 50(10) 

models using extra singlet neutrinos (Ti,i = 1, 2, 3) and the Higgs representations 126h© 
126h as well as 16h © 16h- A singlet neutrino, which we show can be as light as 10^-10^ 
"^ I GeV, decays through its small mixings with right-handed (RH) neutrinos creating a lepton 

asymmetry which is explicitly shown to be flavor dependent. While the doublet vacuum 
expectation value (vev) in 16h triggers the generation of desired mixings, it also induces 
r~| I a large RH triplet vev that breaks the left-right intermediate gauge symmetry and gives 

large right-handed neutrino masses. Manifest unification of gauge couplings and gener- 
^Sl ■ ation of heavy RH neutrino masses are achieved by purely renormalizable interactions. 

^ I The canonical (Type-I) see-saw contributions to the light neutrino mass matrix cancel 

QQ I out while the Type-H see-saw contribution is negligible. Determining the parameters of 

^^ ■ the dominant inverse see-saw formula by using the underlying quark-lepton symmetry 

^*j , and neutrino oscillation data, we show how leptogenesis under the gravitino constraint is 

C^^ I successfully implemented. New formulas for the decay rate and the asymmetry parameter 

^^ ' are derived leading to baryon asymmetry within the observed range without invoking a 

resonant condition on RH neutrinos. The model is found to work for hierarchical as well 
as inverted hierarchical light neutrino masses. Testable predictions of the model are RH 
doubly charged Higgs bosons which may be leptophilic and accessible to the Tevatron, 
LHC or a linear collider. In a model-independent manner, the Drell-Yan pair production 
5^ I cross section at Tevatron or LHC is shown to be bounded between 59%-79% of their left- 

handed counterparts with same mass. In contrast to single-step breaking SUSY GUTs, 
which predict a long proton lifetime for the decay p ^ e'^ir^, here this lifetime is substan- 
tially reduced, bringing it within one order of the current experimental limit. 
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I Introduction 

5*0(10) yy with supersymmetry (SUSY) has been at the centre of attention for a number 
of attractive feaures. It contains just one right-handed (RH) neutrino per generation in its 
spinorial representaion 16. With Pati-Salam [2] and left-right gauge symmetries [3] as its 
subgroups, in addition to unification of the three forces of Nature, it predicts high scale 
unification of quark and lepton masses |1] and has the potential to explain the origin of 
parity (= P) and CP violations. Using the Higgs representations 126h © 126h and IOh, 
it reproduces the small masses and large mixings of neutrinos through Type-I and Type-II 
see-saw mechanisms and their extensions [5], [6l 13 E]- It has been also shown that all the 
fermion masses can be fitted through SUSY 5*0(10) by using suitable Higgs representations |8]. 
Another interesting aspect of the theory is that the observed tiny amount of matter-antimatter 
asymmetry of the universe can be naturally explained through leptogenesis [H] and sphaleron 
effects [101 [m [12]. 

In these theories neutrino masses indicated by oscillation data require the canonical see-saw 
scale of right-handed neutrinos to be in the range of Mr ~ 10^^ - 10^^ GeV. This also sets the 
scale for the masses of associated Higgs triplets carrying B — L = ±2. This scale of neutrino 
mass generation is high in models with variants of the canonical see-saw [71 |8l [12] as well. 

It is well known that the scale of leptogenesis through right-handed neutrino decays and canon- 
ical see-saw is constrained from below leading to the lower bound on the lightest RH neutrino 
mass Mn^ > 10^ GeV [13j. This in turn requires the reheating temperature of the universe after 
inflation to be at least Thh ~ 10^ GeV. On the other hand, big-bang nucleosynthesis in SUSY 
theories sets a severe constraint on the gravitino mass and the reheating temperature leading 
to the upper bound Tjih < 10^ GeV [H]. While thermal leptogenesis in SUSY 5*0(10) with 
high see-saw scale easily satisfies the lower bound, the tension with the gravitino constraint is 
manifest. 

Independent of quark-lepton unified theories, the question of baryogenesis via leptogenesis has 
been addressed in the context of the Standard Model (SM) and the Minimal Supersymmetric 
Standard Model (MSSM) [15] where freedom in the choice of Dirac neutrino Yukawa couplings 
permits fine-tuning them to very small values. In most of these models TeV scale resonant 
leptogenesis [16] is realized by degeneracy between right-handed neutrino masses. A major 
difficulty in having low-scale leptogenesis in SUSY 5*0(10) is the absence of such freedom 
because the underlying quark-lepton symmetry requires these Yukawa couplings to be of the 
same order as the corresponding up-quark Yukawa couplings. This latter difficulty persists even 
in some non-canonical see-saw models and several attempts have been made to bring down the 
scale of leptogenesis [17] . 

Another difficulty in renormalizable SUSY 5*0(10) arises from the gauge coupling unification 
constraint and the need for an SU{2)ji x U{1)b-l breaking intermediate scale that generates 
RH neutrino masses through renormalizable Majorana type interactions. It has been found 
that manifest unification of gauge couplings is spoiled in the presence of Higgs triplets of 
126h © 126h with intermediate symmetries such as SU{2)l x SU{2)r x U{1)b-l x 5f/(3)3c 



or SU{2)l X SU{2)fi x SU{4:)4c even at scales a few orders lower than the GUT-scale unless, 
in the first case, the LR gauge theory and 50(10) are extended to include 5*4 fiavor symmetry 
[IS] or additional light scalar degrees of freedom are introduced at lower scales |T9l 120]. On 
the other hand, there are a number of models with light right-handed gauge bosons |2Tl |22] in 
which in place of the Higgs triplets with B — L = ±2 there are doublets carrying B — L = ±1. 
In contrast to the above scenarios, here we are interested in SUSY SO (10) with both doublet 
and triplet scalars. 

The Higgs triplets in 126h and 126h representations include doubly charged bosons, A^='=. 
Dedicated searches for such doubly charged scalars are being carried out at the Fermilab Teva- 
tron [23]. Both the statistics and the energy reach are expected to be further enhanced at the 
CERN LHC However, the high see-saw scale SUSY 5*0(10) models will offer no prospects for 
these searches as the corresponding masses are large, M^ > 10^^ GeV, while in the class of low 
intermediate scale SUSY 50(10) models where only RH doublets in 16h © 16h are used near 
the TeV scale [21] |22l |2l] no doubly charged Higgs bosons are present. 

In this paper we address the issues of neutrino masses and mixings, low-scale leptogene- 
sis consistent with the gravitino constraint, manifest unification of gauge couplings through 
renormalizable interactions, and testable experimental signatures of the proposed model at 
the Tevatron, LHC or ILC. We construct the desired SUSY 50(10) model including the RH 
triplets in 126h © 126h as well as the RH doublets in 16h © 16h, and three singlet fermions 
(Tj, i = 1,2,3) [25]. We find that a singlet fermion in the mass range Mt = 10^ - 10^ GeV 
can go out of equilibrium to generate lepton asymmetry; its decay is naturally suppressed by 
small mixing with heavy right-handed neutrinos (Ni). The vacuum expectation value of the 
RH-doublet in 16h (or 16h) responsible for this desired small mixing also induces a large vevoi 
the RH-triplets in 126h (or TMh). This breaks SU{2)ji x U{1)b-l ^ U{1)y leading to large 
RH-neutrino masses through renormalizable interactions. We find that although heavy right- 
handed neutrinos are present in the model, the Type-I see-saw contributions to the neutrino 
mass cancel out as has been observed in the context of the Standard Model or its extension 
[26| |27] . The Type-II contribution is also found to be negligible. The dominant contribution 
to light neutrino masses arises through an inverse see-saw which has attracted considerable 
attention over the recent years [221 l2Hl 129] . 

In an earlier work by us and S. K. Majee it was found that gauge coupling unification with 
threshold-like behavior would be possible through the presence of two non-renormalizable dim.- 
5 operators [30]. Here, without using any dim.5 operators, we obtain manifest unification of 
gauge couplings in the renormalizable theory with asymmetric left-right intermediate gauge 
symmetry {g2L 7^ g2R) operative at any scale between 10^ and 10^^ GeV. Further, while lepton 
asymmetry was computed through solutions of Boltzmann equation [30] with an assumption 
about the asymmetry parameter, in this work we derive new analytic formulas for the decay rate 
and the CP- asymmetry parameter and find that they are explicitly fiavor dependent. We then 
show analytically that when the model parameters estimated using the neutrino oscillation data 
are used in our new formula, the model yields desired values of the CP-asymmetry parameter 
leading to the observed baryon to photon density ratio. In addition, we demonstrate that the 
model is consistently successful for both hierarchical as well as invertedly hierarchical light 



neutrino masses. The model leaves its testable signature at the LHC, Tevatron and ILC 
[31] 132] through doubly charged right-handed Higgs scalars A^"*^ in the mass range of 100 GeV 
to a few TeV. Since the decay mode A^^ — )■ VF^iy^ is kinematically forbidden these Higgs 
bosons are leptophilic and predominantly result in like-sign charged bilepton pairs A^"*^ — )■ l^l^. 
The absence of light A^ , A^ states and also the absence of left-handed bilepton pairs in the 
decays would provide signatures specific to this model which are different from other bilepton 
production modes. 

In a model-independent manner without using any structure function data, we show analytically 
that the Drell-Yan hadronic pair production cross section for these RH Higgs bosons is bounded 
between 59% — 79% of that for a left-handed boson of similar mass. 

It is found that, triggered by low mass RH doubly charged Higgs, at the unification scale the 
GUT coupling lies in the strong but perturbative regime and the gauge-boson mediated proton 
decay rate is enhanced. The lifetime Tp{p — > e^Tr*^) is shorter and remains within one order of 
the current experimental limit; this can be reached by the ongoing or planned proton decay 
searches [33] [341 [35] . 

This paper is organized in the following manner. In Sec. II we present the essence of the 
model. Unification of gauge couplings with left-right intermediate symmetry is examined in 
Sec. HI where we also discuss proton lifetime predictions. Derivation of new formulas for the 
singlet-fermion decay rate and the CP-asymmetry parameter are in Sec. IV along with the 
predictions for the baryon asymmetry. In Sec. V we discuss testable predictions of the model 
at the Tevatron, LHC and ILC where we also provide an estimate of the upper and lower 
bounds on the Drell-Yan pair production cross section. A brief summary and conclusions are 
given in Sec. VI. 



II The Model 

In this section, we present the salient features of the model responsible for explaining neutrino 
masses, mixings, and leptogenesis with testable signature at accelerator energies. We consider 
the following pattern of spontaneous symmetry breaking originating from SUSY 5*0(10), 

50(10) ^^ 5f/(2)z. X 5f/(2)« X ?7(1)b-l X 5f/(3)c X D [^2213p] 
^^ 5f/(2)i X 5f/(2)fi X t/(l)B-L X 5f/(3)c [^2213] 
^^ SU(2)l X U{1)y X 5f/(3)c [Gs ' 



(A/z) 



SU{3)c X U{1)q 



The first stage of spontaneous symmetry breaking (SSB) is carried out by assigning GUT scale 
vacuum expectation values to the $54 of 5*0(10) along the direction singlelo under the Pati- 
Salam group SU{2)l x SU{2)ji x SU{A)c = Qps 12] as well as the singlet direction under the 



-•^The role of this vev is discussed in 



left-right gauge group SU{2)l x SU{2)r x U{1)(^b-l) x SU{3)c = ^2213 in the Qps multiplet 
(1, 1, 15) contained in a $210 ^^ SO{10). At this stage D-parity remains intact and the gauge 

couphngs of SU{2)l and SU{2)r are equal, ql = Qr [36]. The second stage of SSB takes place 

(2) 
by assigning vacuum expectation value to the D-Parity odd singlet also contained in <l>2io '^^ 

S'O(IO). By suitable fine tunings of the trilinear couplings beteen 210 and the 126h © 126h 

or 16h©T6h the right handed triplets Ar©Ar C 126h © 126h and the RH doublets 

Xr © Xr C 16h © 16h are made much lighter compared to their left-handed counterparts. 

By adopting higher degree of fine tuning for the RH triplet compared to the RH doublet, the 

components of the RH triplet pairs can be assigned masses between 100 GeV to a few TeV while 

the RH doublet pairs are kept heavier, but sufficiently lighter than the GUT scale. Although we 

do not ascribe any vev directly to the neutral components of the RH-triplets in 126h © 126h, 

we will find that once a vev is assigned to the neutral component of the RH-doublet in 16h, 

the triplet vev is automatically induced. Smaller is the RH-triplet mass fixed by the D-parity 

breaking mechanism, larger is the induced triplet vev. 

The reason behind such ordering of Higgs masses and vevs becomes transparent once we consider 
the Yukawa Lagrangian near the intermediate scale emerging from 50(10), 

Cy = Yi^Li^R<!> + ftPlr^ijRAn + F'^rTxr + ^lT'^T + H.c. (1) 

where '4)l,r are left- (right-) handed lepton doublets and T the three fermion singlet fields, one 
for each generation. The superscript T, of course, denotes transpose. In the (z/, A^, T) basis this 
will lead to a 3 x 3 mass matrijo with vanishing 11, 13, and 31 blocks. 




M, = {p N' T)^\mi, Mn Mx \ \ N' \ . (2) 

MJ ft J \T J ^ 

Here the N — T mixing matrix arises through the vev of the RH-doublet field with Mx = Fv^, 
where v^ = (x?j)? and the RH-Majorana neutrino mass is generated by the induced vev of the 
RH-triplet with M^ = fvR, with vr = {Ar). The vev oi the weak bi-doublet $(2, 2, 0, 1) C IOh 
of 5*0(10) yields the Dirac mass matrix for neutrinos, m^ = Y{^^). 

While implementing leptogenesis in this model through T decays, the out-of equilibrium con- 
dition requires the mixing with RH neutrinos to be small. This will be naturally obtained if 
Mjv > Mx or if 'Ur> v^. 

Assuming Mn ^ Mx ^ fi,mo, which would be highly desirable for the present model, inte- 
grating out the heavy RH-neutrinos leads to the effective Lagrangian 



C^mass) = -{fi-M];M^^Mx)T^T-mDM^'mliy^iy 

-MxM^^ml fu + h.c. (3) 

Interestingly, the block diagonalization of this mass matrix results in a cancellation among 
the Type-I see-saw contributions and the light neutrino mass rriiy is dominated by the inverse 



^Each entry in this mass matrix is a (3x3) block. 



see-saw and one obtains, 



rriu = -mo [M^V(^x) ^1 "^ 



D; 



Mt 



M 



Mn + Mx M^^ M|. 



(4) 
(5) 
(6) 



It will be shown in the next section that the left-handed triplets are near the GUT scale while 
vr ~ 10^" - 10^^ GeV leading to negligible Type-II contribution for light neutrino masses for 
suitable values of the model parameters. 

To see how the induced few is generated, consider the Higgs superpotential near the intermediate 
scale where all GUT-scale masses have decoupled, 



W = Maj^ArAr + M^^xrXr + ^ArXrXr + ^2ArXrXr- 



(7) 



/-i-O ■ 



Using (x5j) = (x?j) = "^x' (^r) ~ (^r) — "^R which requires Ai = A2 = A, the vanishing F-term 
conditions, F^o = F^o = F = F-o = give 

R Ft '^ R '^ R 



''^ = -^M^ Ma,M^^ = 2AX, M^^ = -2\vr. (8) 

The above equations imply that even though no direct vev is ascribed to A^j or A^, a large 
vev is induced once a direct vev is assigned to x%,i the latter being essential to generate the 
desired N — T mixings. With lighter RH-triplet masses Ma — 100 GeV - 1 TeV, it is possible 
to have vr ~ 10^° - 10^^ GeV for v^ = 10^ - 10^ GeV. Since vr ^ v^, the spontaneous breaking 
SU{2)r X U{1)b~l -^ U{1)y takes place at the higher scale generating large RH Majorana 
neutrino masses M^ 3> Mx leading to small Ni — Tj mixings needed to establish the out-of 
equilibrium conditions for leptogenesis. 



We assume the Majorana Yukawa coupling to be diagonal, M^r 

gives Ni — Tj mixing angles, 

Mx.. 
sm $,ij ~ - 



M, 



N, 



disig{MN„MN„Mx,). This 



(9) 



In the present model, the left-handed triplet pair in 126 © 126 acquire mass near the D-parity 
breaking scale Mp ^ Mr. In conventional models even with the left-handed triplet mass 
~ lO^"' - 10^^ GeV, the Type-II see-saw contribution is comparable to the Type-I contribution. 
In this model the Type-II see-saw contribution to the light neutrino mass matrix is 



mil = /A' 



M^Ma' 



(10) 



where Mp is the D-parity violation scale which is also the left-handed triplet mass. Now using 
10^ - 10^ GeV, Vu = 100 GeV, and Mr ^ Mu = lO^^-^ GeV, we obtain. 



Ma = 1 TeV, v^ 



mil = /A'(10-'" - 10-^^) GeV, 



;ii) 



which is at least seven orders of magnitude smaller than the highest value of hierarchical masses 
obtained from the neutrino oscillation data as proposed in [30] . 

Subject to small RG corrections, the underlying quark-lepton unification in SO{10) approxi- 
mates the Dirac neutrino mass matrix with the up-quark mass matrix. The light neutrino mass 
matrix is constructed using the available data on neutrino masses and mixings with a reason- 
able assumption on the leptonic phase of the PMNS matrix. Our strategy is to determine the 
mass eigenvalues and mixings of fermion singlets as well as their mixings with RH neutrinos to 
implement the leptogenesis scenario through their decays as will be discussed in Sec. HVl 

Before addressing the leptogenesis issue we show in the next section that manifest gauge cou- 
pling unification occurs in SUSY 5*0(10) with Ql^ intermediate gauge symmetry. No nonrenor- 
malizable dim. 5 operators are needed to support the unification idea. 



Ill Unification, high Wr mass, proton hfetime 

Manifest unification of gauge couplings converging to a GUT-scale value in SUSY 5*0(10) 
models having left-right intermediate symmetry has been found possible earlier by inclusion of 
additional scalar degrees of freedom beyond those needed for spontaneous symmetry breaking 
[T9| [20} [21] . More recently this method has been evoked to fit masses of all charged fermions 
and for explaining small neutrino masses with ly^-boson mass even at the TeV scale [22]. In 
[30] unification of gauge couplings was accomplished by using threshold-like contributions of 
two nonrenormalizable dim. 5 operators at the GUT scale. Manifest unification has been also 
found to be possible when both the left-right intermediate gauge symmetry and SUSY 5*0(10) 
are extended to contain 5*4 flavor symmetry [18]. The left-right gauge symmetry in that case 
also has unbroken D-parity as well as unbroken R-Parity down to the intermediate scale. In the 
present model there is no flavor symmetry. D-parity is broken at the GUT scale and R-Parity 
is spontaneously broken at a lower scale by the vev of RH doublets in 16h- In addition the 
model has a testable novel feature of accessible doubly charged Higgs scalars. 

In this section we show how manifest unification takes place with the gauge couplings of ^2213 
converging at the GUT-scale without invoking the effect of any nonrenormalizable operators. 
We also show how the proton lifetime for the decay p — )■ e^vr" is brought closer to the current 
experimental limit [31]. 



III.l Gauge coupling unification 

We assume the superpartners of the SM particles to have masses of the order of a TeV. Usinj 
renormalization group equations (RGEs) for the gauge couplings up to one-loop 

dgi / «i ^ 3 



"^"-lle^'*'- '^2' 



where i ranges over the set of gauge couphngs. Below we hst the particles which, with their 
superpartners, contribute to the Oj coefficients in different energy ranges. 

(i) Mz<fi< Msusy: 

Here the particle spectrum is the same as in the non-SUSY SM with three fermion generations, 

41 19 

"-y = TT^' "2L = — —, flac = -7. (13) 

10 D 

(ii) Msusy < I^ < M^: 

In this range, in addition to the MSSM particle spectrum, we have the doubly-charged Higgs 
bosons left as unabsorbed components of RH Higgs triplets and these modify only the ay 
coefficients compared to the MSSM. 

57 
ay = ^, a2L = 1, a^c = -3. (14) 

5 

Because of relatively larger value of ay = ^ {ay = -y for the MSSM), due to the A='=^ near 
the TeV scale, the f/(l)y coupling grows faster, triggering a tendency of unification at substan- 
tially lower scales. This difficulty is bypassed by embedding G213 into the G2213 intermediate 
symmetry. At the boundary point, the U{1)b~l coupling starts from a lower value while the 
SU{2)ji coupling is higher ensuring unification at the GUT scale. The exact unification of 
all four couplings of G2213 is achieved by introducing additional scalar submultiplets such as 
(7^(3, 0, 1) and Cgi^l, 0, 8) at scales M^ and Mc, respectively. It has been noted earlier that such 
states in the adjoint representations of the standard model subgroups with y = could be nat- 
urally light and arise as continuous moduli states of string theory, playing a significant role to 
reconcile the discrepancy between the GUT scale and the string scale [3H]. In our case these 
submultiplets are contained in the 5*0(10) representations 210 and 45, whereas 6*8(1,0,8) is 
also contained in the Higgs representation 54 C 50(10). Alternatively, every pair of triplet a^s 
can be replaced by a fermion triplet which has been noted to play the role of stable dark matter 
[39] if its mass is low. This fermionic state along with others may be present in nonminimal 
5*0(10) representations 



We will show below that one set of solutions of RGEs needs a pair of triplet scalars (rio- = 2) or 
equivalently a fermionic triplet with mass ~ 100 GeV. In that case, only the scalars C*8(l, 0, 8) 
may be treated as naturally light continuous moduli states of string theory, or, purely from 
5*0(10) point of view, the mechanism of Refs. p^ 122] can be utilized to make them light by 
exploiting the generalized superpotential 



Another pertinent question arises if one wishes to use a pair of moduli states cri(3, 1, 0, 1) under 
Glr- How is the lightness of these states ensured in the context of D-Parity breaking at the 
GUT scale leading to lighter components of RH triplets in 126h © 126// and RH doublets in 
16// © 16//. This question is readily answered by examining the part of the superpotential. 



W = Wi + W2 + Wi 



3 



Wi = Mi26l26//126// + Ai26210//126//126//, 

8 



W2 = Mi616^16h + Ai6210^16h16^, 

Ws = M4545^ + A45210h451, . (15) 

Noting that the singlet under Qps contained in 210 h is D-odd, the RH triplets are made light 
when the parameters M^e and A126 < 210^ > are in the same phase. Similarly the RH doublets 
are made lighter than the GUT scale when Mig and Aie < 210// > are in the same phase. Thus, 
it is clear that the same mechanism also renders (Jl{3, 1, 0, 1) C 45// substantially lighter than 
the GUT scale while keeping aR{l,3,0,l) C 45// heavy when M45 and A45 < 210// > are in 
opposite phase. 

Purely from SUSY 5*0 (10) considerations with standard three fermion generations, the method 
of keeping the relevant Higgs scalars substantially lighter than the GUT scale has been discussed 
in Refs. [HI [22] by exploiting the minimization of the generalized superpotential of Ref.|41]. 

(iii) M^< fi<Mc: 

In this range in addition to the contribution of the particles listed above we include those 
from n„ members of Higgs scalar triplets cr/, (3,0, 1) leading to a2L = 1 + ^n^, and ay = ^, 
osc = — 3 as before. 

(iv) Mc < /x < M/j: 

Over and above the contributions mentioned above, here we include the Uc color octets Cs{l, 0, 8) 
resulting in a^c = — 3 + Sric and ay = -y, a2L = 1 + ^.n^, as before. 

(v) Mr<^i< Mu: 

In the presence of Qlr gauge symmetry we have contributions of all the sub- mult iplets discussed 

above. In addition, from the 126// © 126// and 16// © 16// the following submultiplets must now 

be included: 

$(2, 2, 0, l)©Xfl(l, 2, -1, l)©X/j(l, 2, +1, l)©A/j(l, 3, -2, l)©A/j(l, 3, +2, l)©n.aL(3, 1, 0, 1)© 

ncC8(l,l,0,8). 

(vi) /i > Mr with rio- = "^c = 3 we have 

asL = 33/2, a2L = 7, 02/? = 6, aac = 6 . (16) 

With the above particle content and keeping the possibilities of M„,Mc smaller or larger 
than the intermediate scale Mr, allowed solutions are realised with Mr = 10^ - 10^^'^ GeV, 
Mu = 10^^-^^ - 10^^-^ GeV and a^^ ~ 5 - 10. This covers the desired range Mr = 10^^ - 10^^ 
GeV required to implement viable leptogenesis while satisfying the gravitino constraint. 

For a typical example, the evolution of the gauge couplings and unification at the GUT scale 
are shown in Fig. [1] for which we have obtained 

Mr = 10" GeV, Mu = 10^*^ GeV, (17) 

with aQ^ = 5.3 which is well within the perturbative limit. In Fig. [T]the couplings for SU{2)r 
and SU{3)c are found to be almost ovelapping above the scale Mr because of a fortuitous 

9 



ria 


ric 


(GeV) 


Mc 
(GeV) 


2 


2 


100 


7.16 X 10*^ 


3 


3 


3.8 X 10*^ 


8.67 X 10^ 


4 


3 


8.63 X 10^ 


8.67 X 10^ 



Table 1: The number of o"l(3,0, 1) and 6*8(1, 0,8) submultiplets with their respective mass 
scales which lead to solutions with Mjy = 10^^ GeV, M^ = 10^^ GeV, and a^^ = 5.3. 

identity of their respective beta function coefficients and near equality of the boundary values 
at Mji in this example. The change in slopes at M^ and Mc are clearly noticeable. 
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Figure 1: Unification of gauge couplings with left-right symmetry breaking at M^ = 10^^ GeV 
with n„ = nc = 3 (see text). Below Mu = Mp = 10"*^^ GeV, ql 7^ Qr- The SU{3)c coupling 
(short-dashed line) and the SU{2)r coupling (dotted line) appear to merge for fi > Mr because 
fortuitously, in this example, the beta-function coefficients for both are nearly equal in this 
energy range as are the boundary values of the two couplings at Mr. 



In Table [T] we present several choices of n^- and nc and their respective mass scales for which 
the same values of M^ = 10^^ GeV and Mu = 10^^ GeV are obtained as in Fig. [TJ The same 
results follow when a pair of a^s are replaced by a fermion triplet contained in additional 
5*0(10) representations. These fermions while driving type-Ill see-saw for their appropriate 
mass ranges, may also serve as stable dark matter candidates if their mass is low [391 1^ - 



III. 2 Observable gauge boson mediated proton decay 



There are elegant methods and models to suppress Higgsino mediated proton decay or allow 
both types of decays through dim. 5 or dim.6 operators [^ H5| HI]. In most of the single- 



10 



step breaking models, barring a few [S], neglecting threshold effects, the unification scale is 
M^ = 2x 10^^ GeV with a^^ ~ 25 which imply large values of the lifetime for gauge boson medi- 
ated proton decay, e.g. p — )■ e+7r°, for which the current lower bound is [31] {Tp)expt. > 1-01 x 10^^ 
yrs. Extensive estimations of the decay rate have been made in minimal GUTs and their ex- 
tensions with or without SUSY [331 HOI HH SH]. Upto a good approximation, the decay width 
in the present model can be written as 

/ 4 \ 



^(^ -»''"' - 64^il|l-^i-^(^ + « + ^' 



X 



[iAln + Al,){l + K,\r]- (1^ 



In the above formula an is the hadronic matrix element, rup = proton mass = 938.3 MeV, 
/tt =pion decay constant= 139 MeV, and the chiral lagrangian parameters are D = 0.81, F = 
0.47. The short distance renormalization for relevant dim. 6 operators evaluated with super- 
symmetry from Ml! — )■ MsusY and without supersymmetry from Msusy = 1 TeV — )■ Mz in 
the present model with appropriate anomalous dimensions [16] gives Asl — Asr = Asd = 
2.38. The long distance renormalization factor is known to be Al = 1.25. Noting that 
Aji = AlAsd — 2.98, and Fq = 2(1 + iKdP)^ — 7.6, we then express the lifetime as 



V-\p -^ e+7r°) = (1.0 X lO^Vs.) (- 



0.012 GeV^V /2.98V f^l^ 



OiH J \AjiJ \aG 



xl^V.^.41^1 . (19) 



FgJ \1.3x lOi^GeV 

where we have used an = ^^(1 + D + F) ~ 0.012 GeV'^ as per recent lattice estimations 
|47j . In a number of single-step breaking models or other intermediate breaking models with 
M[/ = M^ = 2 X 10^^ GeV and a^^ ^ 25, the one-loop estimation gives large proton lifetime 
Tp{p — )• e + 7r°) ~ 0(10^^) yrs. which is beyond the experimentally accessible limits of ongoing 
and planned proton decay searches for the p — )■ e~^n^ mode. 

In the present model some of our predictions at one- loop level using eq. (IT^ and RGE solutions 
are given in Table [H Although two-loop and threshold corrections are likely to improve these 
results, at one-loop level itself our predictions on the lifetime are substantially less than a large 
number of single step breaking models except few [H] and other intermediate breaking models 
in conventional SUSY 5*0(10) GUTs. Our model predictions are found to remain within one 
order of the current experimental limit and are likely to be accessible to ongoing and planned 
experiments for proton decay searches [5^ [53] . Out of the two observable model predictions, 
namely the low mass RH doubly charged Higgs (discussed later) and proton decay, if any one 
is first observed experimentally, the observation on the other should follow. 
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Mr 

(GeV) 


Mu 
(GeV) 


«c^ 


Tpij) — )■ e+7r°) 
(yrs.) 


IQii 


1.4 X 10^*^ 


5.3 


1.5 X 10=^^ 


10^1 


2 X 10^" 


4.2 


4 X 10^^ 


W 


2 X 10^6 


4.1 


3.8 X 10^4 


10^3 


2 X 10^6 


3.3 


2.5 X 10^4 


10^ 


1.4 X 10^6 


6.2 


8.6 X 10=^^ 



Table 2: Gauge boson mediated decay lifetime for p -^ e+7r° in SUSY 50(10) with (^2213 
intermediate symmetry as described in the text. 

IV Leptogenesis through singlet fermion decay 



IV. 1 Leptogenesis and canonical see-saw 



In the standard formulation of leptogenesis the lightest right-handed neutrino decays into either 
l~(f)'^ and v(j)^ or into their conjugate channels l^(j)~ and vcfP and the desired CP-asymmetry is 
generated by the interference of the tree level amplitude with one-loop amphtudes (vertex and 
self-energy corrections). Denoting the mass eigenvalue of the zth RH neutrino as MjVj, using 
the canonical Type-1 see-saw formula the decay rate of A^i is, 



T, = ^Mm,{Y},Yd)u 



1 



57rf' 



-miM%^ 



(20) 



where rhi is roughly the lightest left-handed Majorana neutrino mass, Vu the vev of the up-type 
Higgs, and Yd the Dirac-type neutrino Yukawa matrix which, upto RG corrections, is the same 
as the up-quark Yukawa matrix. A net lepton asymmetry is generated when the decay process 
goes out of equilibrium at temperature ~ M^Vi satisfying the condition. 



H{T) = IMgl/'^ 



M. 



PI 



(21) 



where H is the Hubble expansion rate. In the normal hierarchical case the generated CP- 
asymmetry is expressed as. 



ei 



3 Mmi Im[{'rrrj~,m£,) 
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IttvI Mjy^ 



m^D^D)ii 



(22) 



The canonical see-saw mechanism gives rise to the lower bound Mjvj > 2.9 x 10 GeV [T3] . 
Since this exceeds the upper bound on the gravitino mass by several orders, the tension between 
standard leptogenesis and the gravitino constraint in SUSY theories is explicit. 

In a large class of solutions of the present model, examples of which are considered in the 
following subsections, the decay of two of the three singlet neutrinos to Z"^"'" and z/0° (and the 
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charge conjugate states) is kineinatically forbidden. Also, the mass of the remaining singlet 
neutrino is ~ 10^ - 10^ GeV, determined by the neutrino oscillation data, and is substantially 
smaller than that of the RH-neutrinos. Thus, its decay can reconcile with the gravitino con- 
straint provided it goes out-of-equilibrium at temperatures ~ Mt^ and it generates the required 
lepton asymmetry. The out-of-equilibrium condition is found to be naturally achieved due to 
the small mixings of Ti with heavy right-handed neutrinos dictated by the model. 



IV. 2 Model parameters for the inverse see-saw 

In this subsection we show how the parameters needed for leptogenesis are obtained in this 
model using the inverse see-saw formula, neutrino oscillation data, up-quark masses and the 
CKM matrix elements. Apart from the RH neutrino mass matrix M^r assumed to be diagonal 
with the largest element of order vr, eqs. (|4]) - ([6]) contain four mass matrices: the light neutrino 
mass matrix 771^, the neutrino Dirac mass matrix m£,, the N — T mixing matrix Mx, and the 
singlet fermion mass matrix fi. 

We construct mu from the mass eigenvalues (^711^2,3) via the PMNS matrix, Upmns, for which 
we use 612 = 32°, 6^23 = 45°, 613 = 7° and take the leptonic Dirac phase Spmns = 1-0 radian. 

m^ = Upjy.jj^g dmg{mi,m2,m3)UpMNS ■ (23) 

The Dirac neutrino mass matrix is fixed by the underlying quark-lepton symmetry of SO{10). 
Neglecting small RG corrections, it is taken to be approximately equal to the up-quark mass 
matrix. Using the PDG values |18] of the CKM matrix elements, its Dirac phase, and the 
running masses of the three up-type quarks, namely, niu = 2 MeV, tjIc = 1.5 GeV, mt = 171 
GeV, we have, 

mo^Mu = V^^^,^ diag(m„, m^, mt)VcKM, (24) 

where we have used the CKM phase 6c km = 1.0 radian, and the quark mixing angles sin6''' — 
0.2243, sin^lg = 0.0413, and sin^^g = 0.0037. 
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The mass matrices Mx and fi are not constrained by experimental data. To minimize the 
unknown parameters we assume a simple form for the N — T mixing matrix Mx, 



where using eqs. ([T]) and 



/ Mxn 

:x= 

V Mx,, 




Mx,, 



e have defined 




Mx,, = Fijv^ 





(25) 



(26) 



With the knowledge of rriu, mo and Mx we then use the inverse see-saw mass formula - eq. 
(jll) - to obtain elements of the matrix /i for hierarchical as well as invertedly hierarchical light 
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neutrino masses. The /i and Mx matrices are used in eq. (jSD to compute the mass eigenvalues 
of the singlet fermions and their mixings by the diagonalization procedure with, 

f, = 5^t/.,T,-, (27) 

i 

where we denote the mass eigenstates by Tj and the corresponding mixing matrix by U . Thus 
the two inputs matrices Mat (chosen diagonal) and Mx - eq. ( I25l) - completely determine the 
singlet neutrino, Tj, masses and mixings consistent with the data on the light neutrino mass 
spectrum, mixing, and grand unification. 



IV. 3 Analytic formulas for decay rate and asymmetry parameter 

The physical processes responsible for leptogenesis are shown in Fig. |2] where crosses denote 
appropriate N — T mixings. The fiavor dependent decay rate for the singlet fermion Tj of mass 
Mj^ through its mixing with Ni - recall eq. (|9]) - can now be expressed as, 

^T, ^ \Uij? sin^ S,jk{in^D^D)kk- (28) 

jk 



ItivI 



T N 




Figure 2: The tree and one-loop level contribution to the decay of Ti that generate the lepton 
asymmetry. 



We find that, depending on the choices of M^ and Mx, there is a wide possibility for the singlet 
neutrino mass eigenvalues. In particular, the model permits a class of solutions where only one 
state has mass above the l(j) threshold (~ 100 GeV) while two others have masses substantially 
below. Denoting this eigenstate as Ti we discuss leptogenesis through its decay in the rest of 
this paper. Because of the simple assumption on the Mx matrix given in eq. ( I25l) the decay 
rate and the asymmetry parameter are reduced to the forms. 



' Stt ' K. 



2 



{\Uu\Ysm'^n{YlYn)n + m2\?sm^^S2{YlYDh3 

+ {\Ul3\?Sm^^23{Y^YD)22\, (29) 
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Hier. 


(lO^GeV) 


Mx,3 
(lO^GeV) 


(lO^GeV) 


(GeV) 


rrii 
(10-2eV) 


f/ii 


t/l2 


t/i3 


(10-^) 


NH 


1.2 


0.18 


18 


1.11x10^ 
92 

7.7 


1.1 
1.4 
5.2 


-0.992 


-0.125 


0.035 


1.11 


IH 


1.2 


60.0 


1.8 


1.83x10^ 
65 
9 


5.3 

5.4 
2.0 


-0.085 


-0.001 


0.996 


1.40 



Table 3: Sample solutions with one singlet, Ti, at the right mass scale for leptogenesis. For 
both normal (NH) and inverted (IH) hierarchies the masses of the singlet neutrinos and the 
light neutrinos are displayed. M^ = diag(2 x 10^, 5 x 10^", 9 x 10^") GeV has been chosen for 
both cases, k, = 0.40 — 0.55 for agreement with the observed baryon asymmetry. 

where Ki, K2 are modified Bessel functions. Even though Y^, is of the same order as the 
up-quark Yukawa matrix, the smallness of Tt^ compared to the Type-I see-saw case - eq. ( 120|) 
- originates from two sources: (i) Allowed values of Mt^ <C M^^ii = 1,2,3), (ii) sin^^j^ <^ 
1 {j,k = 1,2,3). These two features achieve the out-of-equilibrium condition at temperature 
~ Mt-^ satisfying the gravitino constraint. The asymmetry parameter can be expressed as, 

-3Mt P 



87r g' 






( ?7n|)'sin2eii/Mjv3 - ( U12 fsin^^.^/M^, 


MYuy3*i\' 


+ 


{\Uu\ysm'^23/MN, - ( Ui2\ysm'^,2/MN, 


Im[Y2,Y;^' 


+ 


' i\Uu\f sin' U/Mn, - ( Un\f sin' ^23/Mn^ 


Im[YuYl]' 


Un''sm^U{YlYD)n + \Ui2f sin' U{yIYd)33 


+ 1 


Ur3'sin'^23iY^Yn)22 ■ 







Q 

(30) 

In Table [3] we present for two typical solutions the right-handed Majorana neutrino masses, 
matrix elements of Mx, the mixing matrix U , and the T-particle masses when the light neutrino 
masses are hierarchical or invertedly hierarchical, ei is obtained through eq. f pOj) . 



IV. 4 The baryon asymmetry 



For large departure from equilibrium in the Ti decay, the lepton asymmetry per unit entropy 
at temperature T > Mt^ is |49] 

nL_ ^ Ke^grJ^ _ 45 gr^ 

S S TT^ 



-Kei = 4.33 X lO^'^Kei 



(31) 



27r4 g, 

where k is the efficiency factor and gTi = 2 the number of degrees of freedom of Ti . The entropy 
density s = (2/45)(7*7r^T^ where g^, = 106.75, the effective number of relativistic degrees of 
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freedom contributing to entropy in the standard model. Denoting by Nfj the number of Higgs 
doublets {Nh=^ in this model), the baryon to entropy ratio is: 

s 66 + 13Nh s 79 s ^ ^ 

Noting that s = 7.04n^, where n^ is the photon density, the observed baryon asymmetry is, 

riB = — - lO^^Kei. (33) 



n^ 



This is to be compared with [50] : 

(r/ij)expt = (6.15 ± 0.25) x 10^^°. (34) 

We find that for both the cases (NH as well as IH) the predictions are in agreement with the 
observed value. In Table |3] we have exhibited only two out of a large number of allowed solutions 
with efficiency factors k = 0.4 — 0.5. 

For comparison, in Refs. jUJ |5T] maximal efficiency, k ~ 1, has been considered to obtain the 
requisite baryon asymmetry. In Ref. [51] constraints on the Dirac Yukawa coupling of the RH 
neutrino has been examined and it turns out to be small. In our model the effective Yukawa 
coupling of the decaying particle Ti (instead of Ni) to the l(j) pair is essentially modified by the 
product of two extra factors each of which is a mixing substantially smaller than unity. Thus, 
the effective Yukawa coupling of the decaying singlet neutrino remains small. 

The present model permits a variety of solutions with ei ~ 10~^ - 10~® which match the 
observed baryon asymmetry when the efficiency factors n ~ (9(10~^) — 0{1). 



V Right-handed doubly charged Higgs at colliders 

In this section we briefly discuss how this model can be experimentally tested at high energy 
colliders such as the Tevatron, LHC or ILC. The light doubly- charged Higgs boson provides a 
clear scope for this. We relate the production cross section of the A^^ with that of a A^^ of 
the same masqj. This relationship permits the setting of upper and lower bounds on the pair 
production cross sections at the Tevatron or LHC in a model-independent manner. 

As explained in Sec. [TTl the RH-triplets in 126h © 126h carrying B — L = ±2 are made light 
in this model through the D-parity breaking mechanism with the component masses M^,^ — 
100 GeV to a few TeV. This enhances the induced vev, vr, resulting in high scale LR gauge 
symmetry breaking with large Wj^ , Zji gauge boson masses and heavy RH Majorana neutrinos. 
The largeness of the RH Majorana masses in the model lead to naturally small T-N mixings 
essential to satisfy the out-of-equilibrium condition for the T-decay to drive leptogenesis. 



^Iii our model the A^ is very heavy due to D-Parity breaking. 
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Because of large masses of the RH gauge bosons, myy±,mzji ~ ^R ~ ^^ ' ^^ GeV, the decays 
A^^ — )■ W^^W^^ are kinematically forbidden and so only leptonic decays are possible. The 
5*0(10) invariant Yukawa interaction with fermions given in eq. ([1]), fij.lQi.16j.12QH, makes the 
doubly charged Higgs boson leptophihc, its only decay modes being A^^ — )■ l%lR.{i,j = 1, 2, 3). 

Due to their heaviness, Wji and Zji also play no role in A^^ production at the Tevatron or 
LHC The production will dominantly be through the electromagnetic or Z*-exchange Drell-Yan 
mechanism. 

The couplings fij determine the branching ratios of the A^ to the different leptonic final 
states. Although for the sake of economy and simplicity we have chosen fij to be diagonal in 
the previous section, our choice has been guided by negligibly small values of the non-diagonal 
elements suggested by current limits on lepton flavor violating decays such as yU — t- 3e and 
r — )■ 3e. Out of a large number of possible solutions, the two sets given in Table |3] have 
Mn^ = 2 X 10^ GeV, M^^ = 5 x 10^° GeV, Mn^ = 9 x 10^° GeV, which for vr ~ 10^^ 
GeV corresponds to /i = fee = 0.0002, /2 = /^^ = 0.5, /s = frr = 0.9. Keeping in mind the 
wide classes of solutions permitted in this model, we will discuss possible implications for /i = 
0.0002 - 0.001, and /a < /a. 

The mass ordering of the W^, Zr, and A^^ discussed above is specific to this model. Interest- 
ingly, low mass doubly charged Higgs bosons with similar interactions have been shown to be 
generic in a class of SUSYLR models with Mwji > 10^ GeV which require non-renormalizable 
terms in the superpotential [521 |53l |5l]. Within the non-SUSY left-right model prospects of 
Drell-Yan pair production of LH doubly charged Higgs bosons at high energy colliders and their 
detection [53] and the impact of QCD corrections thereon [SB] have also been investigated. 



V.l Bounds from muonium-antimuoniuni conversion 

Muonium(M)-Anti-muonium(M) conversion, /x+e^ — ;■ /x^e+ can be mediated by A^^ giving 
rise to an effective coupling [57] . 

Experimental searches for this transition yields the upper bound [58] 

Gm-m < 3 X IQ-^Gf, (36) 

where Gp = 1.17 x 10"^ GeV~^ is the Fermi coupling. Combining eqs. fl35l) and fl36l) gives 

Ma > ( ^^^^^ — 1 . (37) 

V12 X 10-3v^G^; ^ ^ 

For example, choosing /i ~ /2 ~ 1, a condition applicable to quasi-degenerate RH neutrinos 
but which is outside the presently allowed solutions, we obtain Ma > 2.24 TeV which is beyond 
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the Tevatron limit, but within the LHC range. But when we use the class of solutions which 
permit /i ~ 0.0002 and /2 — 0.5, we obtain Ma > 22.4 GeV. This is not inconsistent with the 
experimental search limit reached by DO and CDF collaborations at Fermilab Tevatron with 
Ma > 112 GeV/c^ and Ma > 127 GeV/c^, respectively [Ml EO]. 



V.2 Drell-Yan pair production at LHC or Tevatron 

Previous searches at LEP have already excluded Aj^ below 97 GeV/c^ [6T]. At hadron colliders 
the doubly-charged Higgs boson will be dominantly created through pair production via the 
basic Drell-Yan process qq — )■ 'j*,Z* — )■ A^"''A^^ ,{i = L,R). At the quark level, the cross 
section depends only on the quantum numbers and mass of the doubly charged scalars. In the 
present model with supersymmetry and purely renormalizable interactions, the only allowed 
decay mode is. 



A^ 



~^ ''R''R 



{31 



and its conjugate. At the Fermilab Tevatron or the LHC, the production of the doubly charged 
boson will be through 



pp(p) ^ (7*, Z*)X ^ A++A^-X -> Ipp'nl'nX 



(39) 



where for the dominant modes in our model 1,1 = r, fi since fi <^ f2 — fs — 1. 

pair production cross section through 7* and Z* exchange is expressed 



The parton level A," 
as 



CTi 



QlQl + 



's-MlY + TlMl 






L,R. 



(40) 



Here s = Q^ = ts is the square of the cm. energy of the colliding quark-antiquark pair and 
r the product of momentum fractions carried by them, a is the fine-structure constant at the 

relevant energy scale and /3j = w(l — 4M^ys) is the velocity of the doubly charged Higgs, 

A^^, produced in the colhsion (i = L,R). Qq^ = {Isq — 2(5gSin^6'H/)/(2sin6'iy cos6'vi/), (7gA = 
Isq / {2 sin 6 w cos 6 w) for the quark g and ^fA^i; = {hA,~QA,s'^)/{2swC^a). Qq{QA,) is the electric 
charge number of the quark q (Higgs Aj) and hq^hAi) the third component of SU{2)l isospin 
for q (Aj). 



It is clear from eq. (HOj) that for 



s > 4Ml > M|, 



(41) 



the S^ become independent of s and the momentum fractions carried by the quarks leading to 
S? ^ S^ = [QlQl^ + {g'^^ + gl)gl, + 2Q,Q^^g,,g^^,] .{i = L,R). (42) 

Noting that Jsa^ = 1(0) for i = L{R) we obtain from eq. ( l42l) 

S^ = 0.59S^, % = 0.79SJ. (43) 

These relations translate to upper and lower bounds on the production cross section of the 
A^^ which do not depend on the proton structure functions or the model-origin of these Higgs 
bosons as long as the gauge symmetry at the electroweak scale is the standard model. Thus, 

0.59aL <aR< 0.79(Tl. (44) 

The above bounds for Drell-Yan pair production relate the cross section for RH doubly charged 
Higgs at LHC or Tevatron with that for their LH counterparts with the same mass provided 
the mass is > 150 GeV. 

Using these relations and the published results for A^^ production one can estimate the number 
of events in this model. For examplcl, at the LHC, for m^±± = 200 GeV (1 TeV) the Drell-Yan 
production cross section in fb is 49.4, 96.4, 169.5 (0.004, 0^04, 0.14) for ^/s = 7, 10, or 14 TeV, 
respectively. Using eq. ( H^ and assuming an integrated luminosity of 30 fb~^ for m^±± = 200 
GeV one would expect 1022, 1995, 3509 events for the three cases. If m^±± = 1 TeV then for 
^/s = 10 TeV one would require about 200 fb~^ integrated luminosity for a 5-event signal. 



At the Tevatron with ^/s = 2 TeV, a^j; ^ 12 — 16 fb for M^^ = 150 GeV and with an integrated 
luminosity of 350 pb^^ the predicted number of events is nearly 4 — 6. With an acquired 
integrated luminosity of 10 fb^^, the mass reach of up to M^,^ = 300 GeV can be achieved 
where a/? ~ (5 — 7) x 10"^ fb. 

In our model since /i ^ /2 ~ /s , the dominant decay modes of the produced pair would be 
through the following four lepton channels every one of which would be almost equally likely 

^R ^R ~^ '^R '^R '^R '^R ' f^Rf^Rf^Rf^Ry '^R ^R f^Rf^Rj '^R '^R f'Rf^R- 

The standard model backgrounds for such production processes have been discussed in [55] . 
The signal event should have negligible missing pt- Moreover, the two pairs of like-charged 
leptons are constrained to each have the invariant mass equal to ttia^- These criteria and a 
judicious cut on the /"'"/" pair invariant mass to remove ZZ contributions effectively removes 
the entire background. 

Unlike the Tevatron and the LHC where the doubly charged bosons are pair-produced, at 
proposed muon colliders resonant production of these bosons could take place if fi~fi~ colliders 
are arranged. The singly produced A would decay via 2/i or 2r channels providing the 
cleanest signals for these bosons. In contrast to a large class of asymmetric left-right models 
where the decay of the right-handed doubly charged bosons could proceed via kinematically 
allowed channels such as A^"*" — )■ W^/j^iy^, W^A~I^, A^A^, this model allows decay only in the 
bilepton channel providing a signature of its genuine leptophilic property. 



We thank Anindya Datta for providing these numbers. 

19 



VI Summary and conclusions 

In summary, we have implemented flavor- dependent leptogenesis through the decay of singlet 
fermions with masses M^^ = 10^ - 10^ GeV in SUSY SO{10) while satisfying the gravitino 
constraint. The left-right intermediate symmetry is at a high scale corresponding to Wpt and 
Zr masses larger than 10^^ GeV. This has been made possible by using the RH triplets in 
126h © 126h as well as the RH doublets in 16h © 16h- Not only is the singlet fermion-RH- 
neutrino mixing generated by the vev of 16h but also this mixing becomes naturally small 
through the large vev of the RH-triplet induced by the doublet vev. In addition to obtaining 
renormalizable mass for the RH neutrino through this mechanism, manifest unification of gauge 
couplings is also achieved purely by renormalizable interactions and, thus, non-renormalizable 
dim. 5 operators used earlier are dispensed with. In contrast to the earlier attempts where 
an assumed value of the CP-asymmetry parameter was shown to yield the lepton asymmetry 
numerically, in this work we have derived and suggested new analytic formulas leading to the 
correct asymmetry parameter and the observed baryon to photon density ratio. We have found 
that both the decay rate and the CP-asymmetry are explicitly flavor dependent. Whereas our 
previous work [30j required a normal hierarchy, here we have also found successful implementa- 
tion in the case of inverted hierarchical neutrino masses. Unlike a host of low-scale leptogenesis 
models, this model works with hierarchical heavy RH-neutrino masses and no resonant con- 
dition with extreme degeneracy among them is needed. A decisive test of the present model 
would be through the detection of doubly charged Higgs bosons A^^ — )■ /^/^ at the Tevatron, 
LHC, or a future muon collider. The model provides an example of truly leptophilic doubly 
charged Higgs bosons. Without using parton-density distribution functions we have also shown 
in a model independent manner that the pair production cross sections for RH doubly charged 
Higgs at Tevatron or LHC energies are bounded between 59% — 79% of their LH counterparts 
with same masses. 

As the unification, triggerred by low-mass doubly charged bosons, occurs with a large (but 
perturbative) unified gauge coupling, the decay hfetime Tp{p — )■ e+7r°) is substantially reduced 
compared to conventional SUSY GUTs and remains within one order of the current experimen- 
tal limit. This is likely to be accessible to the ongoing and planned proton deacy searches. The 
model appears to be rich in dark matter candidates which will be investigated elsewhere. 
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